Disorder effects on the spin-Hall current in a diffusive Rashba two-dimensional 

heavy-hole system 



S. Y. LivQ and X. L. Lei 

Department of Physics, Shanghai Jiaotong University, 1954 Huashan Road, Shanghai 200030, China 

(Dated: February 2, 2008) 

We investigate the spin-Hall effect in a two-dimensional heavy-hole system with Rashba spin-orbit 
coupling using a nonequilibrium Green's function approach. Both the short- and long-range disorder 
scatterings are considered in the self-consistent Born approximation. We find that, in the case of 
long-range collisions, the disorder-mediated process leads to an enhancement of the spin-Hall current 
at high heavy-hole density, whereas for short-range scatterings it gives a vanishing contribution. This 
result suggests that the recently observed spin-Hall effect in experiment is a result of the sum of the 
intrinsic and disorder-mediated contributions. We have also calculated the temperature dependence 
of spin-Hall conductivity, which reveals a decrease with increasing the temperature. 
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I. INTRODUCTION 

Spin-Hall effect has been extensively studied in semi- 
conductors with spin-orbit (SO) coupling due to its po- 
tential applications in the emerging field of spintronics. 
It refers to an appearance of the net spin current flow- 
ing perpendicular to the applied dc field. This trans- 
port phenomenon was studied by Dyakonov and Perel 
in 197li and by Hirsch more recently^ Therein, the au- 
thors demonstrated that the spin-dependent scattering 
between electrons and impurities causes the electrons 
with opposite spins to veer to different sides of the con- 
sidered systems. Obviously, such spin-Hall effect strongly 
relies on the electron-disorder collision and hence, lately, 
has been named as extrinsic spin-Hall effect £ 

Recently, an impurity- free spin-Hall effect, namely in- 
trinsic spin-Hall effect, has been proposed by Sinova et 
al. in two-dimensional (2D) electron systems^ and Mu- 
rakami et al. in p-doped bulk semiconductors^ respec- 
tively. This spin-Hall effect is associated with the dc- 
field-induced transitions between the spin-orbit-couplcd 
bands and contributes from all occupied electron states 
below the Fermi energy^ In 2D electron semiconductors 
with Rashba and/or Dresselhaus spin-orbit coupling, the 
value of intrinsic spin-Hall conductivity has been found 
to be a universal value e/8wm& Further, in 2D Rashba 
heavy-hole (HH) systems, Schliemann and Loss reported 
that the intrinsic spin-Hall conductivity starts out at 
value 9e/87r and increases with increasing the Rashba 
coupling. 7 

However, the disorder can strongly af- 
fect the spin-Hall effect, especially in 2D 
semiconductors! 5 ! 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 : 18 There exists 
another collision-related spin-Hall current connecting 
to the electron states near the Fermi surface. 5 It is 
disorder-mediated but independent of the impurity 
density. In diffusive Rashba 2D electron systems, this 
disorder-mediated spin-Hall effect leads to a complete 
cancellation of the to tal spi n-Hall current for both the 
short- (Refs. Illlfl2lll3lll4ll5j) and long-range disorders^ 
The same conclusion has been obtained even in the 



presence of both the linear in momentum Rashba and 
Dresselhaus SO couplings ^LiLi2i 

Experimentally, observations of spin-Hall effect have 
been reported recently by Kato et ali& and by Wundcr- 
lich et alm& The latter work is of special interest, since 
there the spin-Hall effect has been found in a relatively 
clean 2D heavy-hole sample with Rashba spin-orbit cou- 
pling. Obviously, in such a system the observed spin-Hall 
effect cannot be understood as an extrinsic one. To in- 
terpret this experimental finding, it immediately arises 
the question whether the sum of the intrinsic and the 
disorder-mediated contributions to spin-Hall current also 
vanishes in the Rashba 2D heavy- hole semiconductors. 
By means of Kubo formalism, Bernevig and Zhang found 
that the disorder has no effect on the spin-Hall effect and 
the observed spin-Hall effect is a pure intrinsic one^i 

Based on a nonequilibrium kinetic equation approach, 
we demonstrate in this paper that the vanishing of the 
disorder effect on spin-Hall current in 2D HH systems oc- 
curs only for short-range hole-impurity scattering. When 
the disorder becomes long-ranged, the contribution from 
disorder-mediated process to the spin-Hall current can 
not be ignored. It has the same sign as the intrinsic one at 
high HH density, but opposite for low density. The com- 
parison between theory and experiment indicates that 
the observed spin-Hall effect is the result of the sum of 
the disorder-mediated and intrinsic contributions. We 
also find that the spin-Hall conductivity decreases with 
the rise of the temperature. 

The paper is organized as follows. In Sec. II the nonin- 
teracting Green's functions in spin basis and the kinetic 
equations for nonequilibrium distribution functions are 
presented. In Sec. Ill we analyse the disorder effects 
on spin-Hall conductivity for short- and long-range hole- 
impurity scatterings, respectively. Finally, we conclude 
our results in Sec. IV. 
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II. FORMALISM 

In bulk semiconductors like GaAs, the heavy- and 
light-hole (LH) bands with total angular momentum 
j = 3/2 are degenerate at the band edge. However, 
in quasi-two-dimensional semiconductors, the additional 
confinement along the growth direction of heterostruc- 
ture yields a splitting between these bands and the de- 
generacy is lifted. As a result, the HH and LH subbands 
are separated. 

We consider a two-dimensional HH system, where only 
the lowest heavy-hole subband is occupied. Note that 
this condition can be satisfied when the quasi-2D system 
is sufficiently narrow and the density and the tempera- 
ture are not too high. In this way, the effective Hamilto- 
nian for a single heavy-hole subjected a spin-orbit inter- 
action due to the structural inversion asymmetry, can be 
written ao 7 i 22 i 23 i 24 i 25 i 26 i 27 



(1) 



with 2D momentum p = (j>cos(/> p ,psin</> p ), fr± = fr x ± 
a y , and p± = p x ± ip y . fr = (fr x ,fry,frz) are the Pauli 
matrices and m is the hole effective mass. Note that 
this Hamiltonian is obtained through the decomposition 
of the HH-LH coupling by third-order Lowdin perturba- 
tion theoryJ^ The parameter a is proportional to the 
Rashba spin-orbit coupling constant, and relies on the 
separation between HH and LH subbands and the Lut- 
tinger parametersi 2 ^ 2 ^ We can see from Eq. that the 
spin-orbit interaction is cubic in momentum. There ex- 
ists another type of the cubic SO coupling, namely the 
Dresselhaus spin-orbit coupling in 2D electron systems. 
The spin-Hall effect for such SO interaction has already 
been studied in Ref.lla. 

It follows from Eq. that the Green's function of non- 
interacting particles has the form 

G a (u;,p) = n«GSf(w, P ) + n( 2 >G£ 2 >, P ), (2) 



with Gq^ 



z tl {p)±i5]- 1 (M=l,2). Here, 

„2 



EL 

2m 



n 



(-l)A 



(p 3 -(fr 



+ ) a f3 -pAa-) af3 j 



(3) 



, (4) 



and Ehh = otp 3 . The operator fl^ represents the projec- 
tion onto the states with a definite helicity. 

By means of a local unitary transformation 



^ (P) ~ L/2 ( *e 3 ^ 



(5) 



the Hamiltonian can be diagonalized as H = 
U + HU — diag(ei(p), e 2 (p)), and hence the noninter- 
acting ^Green's function also becomes diagonal Gq" = 
U + G U = diag(Gg 1 a , Go' 2 a ). This implies that the spin- 
orbit interaction can lead to further lifting the degener- 
acy of the HH band. In result, two spin-orbit-couplcd 
bands with dispersion relations £ M (p) are formed. Note 
that above transformation corresponds to a change of the 
basis of the Hilbert space from a spin to a helicity one. 
In the following study, the spin-Hall effect will be carried 
out in the helicity basis. 

We assume a weak dc field E applied to the 2D HH 
systems along the x direction. We are interested in the 
z-direction-polarized spin current flowing along the y di- 
rection, i. e. the spin-Hall current. It is known that the 
single-particle spin current operator can be defined by 7 - 



3 1 



(6) 



where the factor | reflects the angular momentum of the 
heavy-hole and the j y is the electric current operator. 
By taking the unitary transformation, j* becomes non- 
diagonal and the net observed spin-Hall current can be 
expressed as 



J y=J2^P^ (p) + &i(p) 
p 



3 P y 



(7) 



The p a /3 (a, j3 = 1, 2) are the elements of the distribution 
function p(p): p a p{p) = -i J ^G^(p,uj) with G< be- 
ing the less Green's function. The spin-Hall conductivity 
is given by a sH = JfjE. 

To carry out the spin-Hall conductivity, it is necessary 
to derive the kinetic equation for less Green's function. 
Following the procedures described in Ref .0, we find that 
p obeys the equation 



3ieE 

eE • V p /5(p) + — — • V p 4> P [p, fr x ] + i£HH[f>, = 



dl 
dT 



(8) 



At the same time, the self-energies for the hole-impurity interaction with isotropic scattering matrix V(|p — k|) can 
be written as 



£ r ><(p,uj) = in i 5^|V(|p-k|)| 2 {ai6 r ' < +a 2 fr x G r ><fr x +ia 3 {fr x ,G r ><}} 



(9) 
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with rii being the impurity density and cii(i = 1, 2, 3) the 
factors associated with the directions of the momenta, 
1 + cos(30 p - 30k), 02 = 1 — cos(30 p - 30k), 03 = 



a 1 



sin(30 p — 30k)- The relaxation term 

scat! 

general form^i 

Similar to the case of 2D electron systems, the lin- 
earized kinetic equation can be further simplified and its 

I 



solution is found to comprise two terms. The first one is 



takes a Pu Cp) =P2i (p) 



3eE ■ V 



^-^ {n F [ei(p)] - n F [e 2 (p)}} 
( 10 ) 

and the real part of the second one takes the form 
(2) 

R" e Pi2 (p) — C(p) e £- sm 0p w ith function £(p) 



C(p) 



4ap 3 



E 



(-1)" dn F (E) 



T4fj.fj. 



dE 



*/.Mp)]- 



(-l)^ 1 <9n F (£) 



T4 W 



as 



*p[ e /»(P)] 



(11) 



The functions <&^(E) connect to the diagonal distribution 
functions 



dn F {E) 



dE 



$ M [e„(p)]e£cos0 p (12) 



and can be given by the coupled equations, 
cMp) _ $^M;p)] , $„[e M (p)] $ mMp)] 



<9p 



T 2>. 



, (13) 



with /i = 3 — (X. In these equations, denote the 
different relaxation times 

— = 2™, \ V (\P - k|)| 2 A^(0 k - P ,P, k), (14) 



where the 
h\nu{<t>>P,k) 



functions A^„(0,p, k) are defined as 
= |[1 + cos(30)](l - cos0)<5(e MP - e„ fc ), 
= |[l-cos(30)]5(£ w -e 1/fe ), A 3lll/ ((j>,p,k) = 



icos0[l — cos(30)](5(£ p 
i sin0sin(30)5(e 



, fe ) and A 4M „(0,p, ft) 



The first term in the solution of the linearized kinetic 
equation gives rise to an intrinsic spin-Hall conductivity 



'sfl 



r (1 > 



-9e 



lQirma 



^{n F [ei(p) 



(A - n F [e 2 {p) - A*]}, 
(15) 

in agreement with the previous studyi This part of spin- 
Hall conductivity comes from the interband transition 
processes in which all holes below Fermi surface join. It is 
associated with the energy separation between two spin- 
orbit-coupled bands in helicity basis and reveals an in- 
trinsic character: it is independent of any hole-impurity 
scattering. Physically, the dc field can cause an elastic 
transition of a hole from one band to another one when 
this hole gains an additional momentum from the exter- 
nal dc field. 

The second term in the solution is associated with the 
transport process and only the hole states near the Fermi 
surface contribute. It relates to the disorder collision in 



a surprising way: its contribution to the spin-Hall con- 

(2) 

ductivity a sI ^ is independent of the impurity density rii 
but depends on the above-defined relaxation times and 
hence the form of scattering matrix. We can understand 
the origin of the er^ as follows. When an external dc field 
is applied, the holes should participate in the longitudi- 
nal transport, leading to a diagonal distribution propor- 
tional to the inverse of the impurity density. At the same 
time, these perturbative holes also experience the impu- 
rity scattering, yielding the interband polarization. In 
result, the nondiagonal distribution becomes impurity- 
density- independent. Here the impurity plays only an 
intermediate role. 

We should note that this disorder-mediated mecha- 
nism is physically identical with the side-jump mecha- 
nism studied recently in the investigation on spin-Hall 
effect m& Contributions to spin-Hall conductivity from 
both mechanisms are independent of the impurity den- 
sity but collision-related. At the same time, they all 
are associated with the hole (or electron) states near the 
Fermi surface. However, formally, these two mechanisms 
look different because of two different approaches. The 
side-jump process is carried out for a spin-orbit inter- 
action involved in the electron-impurity scattering. It 
corresponds to a lateral displacement of the center of the 
wave-packet during the scattering and hence connects to 
the scattering-dependent term of the current operator. 
However, in our study, since the spin-orbit coupling is 
included in the free-hole (or free-electron) Hamiltonian, 
the current operator is independent of the hole-impurity 
collision. 

The fact that the spin-Hall current consists of two 
parts, is similar to the well-known result of Stfeda^ 
as well as the recently obtained conclusion in the con- 
text of anomalous Hall effect (AHE). 30 In the 2D elec- 
tron systems magnetically or with magnetization, the off- 
diagonal conductivity usually comes from two terms, one 
of which is due to the electron states near the Fermi en- 
ergy and the another one is related to the contribution 
of all occupied electron states below the Fermi energy. 
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III. RESULTS AND DISCUSSIONS 

First, we consider the spin-Hall effect for the short- 
range disorder. In this case, the holc-impurity colli- 
sion matrix has a simple momentum-independent form 
V(|p — k|) = u and the scattering is described by a single 
relaxation time t = 1/riiU 2 . It can be seen from Eq. 
that the disorder-mediated contribution to spin-Hall cur- 
rent relies on function A^ pv proportional to sin <f> sin (30). 
When the angle integration in Eq. JJJ) is performed, the 

(2) 

vanishing a „ is obtained. 

I — I 

This result agrees with that of Rcf. 21, in which 
the vanishing of disorder-related spin-Hall conductivity 
comes from an analogous angle-integration. However, for 
2D electrons with Rashba spin-orbit coupling, the func- 
tion depends on sin 2 cf>, leading to a nonzero er^- 5 

For long-range disorders, the additional momentum- 
dependence of scattering matrix produces rich novel phe- 
nomena. We have performed a numerical calculation to 
investigate the long-range disorder effect on the spin-Hall 
conductivity in a 2D GaAs/AlGaAs based heavy-hole 
system. Consider a Coulomb interaction between the 2D 
heavy-holes and the charged impurities located at a dis- 
tance s = 500 A from the 2D plane: V(p) ~ e~ sp I(p)^ 
I(p) is the form factor. In calculation, we take the ef- 
fective mass m — 0.27m e (Ref. |20[) and the coupling 
constant ma = 5 A. 

In Fig. 1 the total and disorder-mediated spin-Hall con- 
ductivities are plotted as functions of heavy-hole density 
at zero temperature. The intrinsic spin-Hall conductivity 
is almost independent of HH density and has the value 
9e/87r, in agreement with the result in Ref.0. It is evident 
from Fig. 1 that, unlike the case of short-range disorder, 
here the a sI ^ and <j s h decrease with descending HH den- 
sity, even becomes negative for n p < 2.5 x 10 10 cm~ 2 
due to the sign change of the quantities . 

For HH density of order of 10 10 cm~ 2 , the Fermi 
energy and the temperature (less than 4.2 K) become 
comparable. Hence, the HH spin-Hall effect should be 
strongly influenced by temperature. The calculation in- 
dicates that with ascending the temperature, the intrinsic 
spin-Hall conductivity begins with the zero-temperature 
value 9e/87r and descends for heavy-hole densities n p in 
the range 1 — 100 x 10 10 cm~ 2 . However, as shown in 
Fig. 2(a), for the disorder-mediated spin-Hall conductiv- 
ity, there exists a critical HH density n pc about equal to 

4 x 10 10 cm~ 2 . With a rise of the temperature, the cr^ 
increases for n p < n pc , while it decreases for the opposite 
case. However, one can see from Fig. 2(b) that the total 
a S R always decreases with increasing temperature. This 
implies that the descending intrinsic spin-Hall conductiv- 
ity plays a dominant role at low HH density. 

Experimental study by Wunderlich et al. has been per- 
formed for heavy-hole density n p = 2 x 10 12 cm~ 2 iSI At 
such high HH density, the value of spin-Hall conductiv- 
ity remains almost unchanged for temperatures less than 
4.2 K. At the same time, the disorder-mediated process 
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FIG. 1: Density dependence of disorder-mediated and total 
spin-Hall conductivity in a Rashba 2D heavy-hole system. 
The lattice temperature T — OK and the coupling constant 
ma = 5 A. 
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FIG. 2: Temperature dependence of disorder-mediated and 
total spin-Hall conductivity in a Rashba 2D heavy-hole sys- 
tem. From bottom to top, the corresponding HH densities 
are n p — 1, 2, 3, 4, 5, 6, and 7 x 10 10 cm -2 , respectively. The 
coupling constant ma = 5 A. 
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plays a positive role and its contribution to spin-Hall con- 
ductivity adds to the intrinsic one. Hence, the 2D HH 
systems with such a high HH density serve as good can- 
didates for the observation of the spin-Hall effect. 

To compare with the experiment, we have calculated 
the spin-Hall conductivity for the HH system with hole 
density equal to the experimental one. The Fermi energy 
is approximately chosen to be 17.5 meV and hence the 
coupling constant is ma = 2.5A. AtT = 4.2K, we obtain 
the total spin-Hall conductivity er s # = 18.2e/87r. This 
result is independent of the impurity density, but relies 
on the forms of the scattering matrix. Although the spin- 
Hall conductivity is not a quantity directly observed in 
experiments, our investigation suggests that the observed 
spin-Hall effect may not only be a pure intrinsic one, as 
claimed in Ref. l2ll rather comes from both the intrinsic 
and the disorder-mediated processes. 

IV. CONCLUSION 

The spin-Hall effect of a 2D heavy-hole system with 
Rashba spin-orbit interaction has been investigated by 
means of a nonequilibrium Green's functions approach. 
We have studied the intrinsic and disorder-mediated 
contributions to spin-Hall conductivity considering both 
the short- and long-range hole-disorder scatterings. For 



short-range collisions, the disorder-mediated contribu- 
tion vanishes. When impurity scattering becomes long- 
ranged, however, the disorder-mediated spin-Hall con- 
ductivity becomes finite and can change the sign. It is 
negative for low HH density and positive for high den- 
sity. The temperature dependence of disorder-mediated 
and total spin-Hall conductivity has also been analyzed. 
The total spin-Hall conductivity always descends with 
rising temperature at different HH densities, whereas the 
behaviors of the disorder-mediated one versus temper- 
ature become hole-density-related. A comparison with 
recent experiment indicates that the observed spin-Hall 
effect is probably a result of the contributions from both 
the intrinsic and disorder-mediated processes. 
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